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Abstract. A new computational scheme for the nonlinear cosmological matter power 
spectrum (PS) is presented. Our method is based on evolution equations in time, 
which can be cast in a form extremely convenient for fast numerical evaluations. 
A nonlinear PS is obtained in a time comparable to that needed for a simple 1- 
loop computation, and the numerical implementation is very simple. Our results 
agree with N-body simulations at the percent level in the BAO range of scales, 
and at the few-percent level up to /c ~ Ih/Mpc at z ~ 0.5, thereby opening the 
possibility of applying this tool to scales interesting for weak lensing. We clarify the 
approximations inherent to this approach as well as its relations to previous ones, such 
as the Time Rcnormalization Group, and the miilti-point propagator expansion. We 
discuss possible lines of improvements of the method and its intrinsic limitations by 
multi streaming at small scales and low redshifts. 
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1. Introduction 

Cosmological perturbation theory (PT) as a tool to study the Large Scale Structure 
of the Universe (LSS) (for a review see [1]) has received considerable interest in the 
recent past. The main motivation is the study of Baryonic Acoustic Oscillations (BAO) 
imprinted in the matter power spectrum (PS), which are one of the main observables 
of present and future galaxy surveys [2l|3llll[5l|6l[71[8l|9l [10]. The goal of these 
measurements is to derive the acoustic scale at the percent level accuracy, in order 
to provide constraints on the Dark Energy equation of state competitive with those 
obtained from measurements of the magnitude-redshift relation of type la Supernovae 
[TT| [T2] . Of course, on the theory side, the matter PS in the BAO range of scales 
(~ 100 Mpc/h) must be computed at the same level. For a ACDM cosmology, the goal 
can be accomplished by means of high accuracy N-body simulations, once the various 
issues related to precise initial conditions, very large simulation volumes, mass resolution 
and time stepping are carefully addressed, as discussed, e.g., in [13]. However, the long 
computational times required by these simulations make it impossible to implement grid 
based or Markov Chain Monte Carlo (MCMC) estimations of cosmological parameters, 
which typically require the evaluation of thousands of PS's. Moreover, if departures from 
ACDM cosmologies are taken into account, for instance by the inclusion of massive 
neutrinos, non gaussian initial conditions, or f{R) theories of gravity, the N-body 
approach is still far from being firmly established at the percent level in the relevant 
range of scales. 

The usage of PT can potentially help in most of these respects. First, being 
based on analytic, or semi-analytic, techniques, computational times are in general 
greatly reduced with respect to N-body simulations. Moreover, although the basic 
formalism is derived for an Einstein-de Sitter cosmology, its extension to ACDM is 
straightforward, and also its formulations in non-standard cosmologies are feasible and 
under control [m [151 [El HZl [13 US]- However, the real boost to these methods was 
given by the papers of Crocce and Scoccimarro (CS) [201 121], who showed how the 
poorly behaving PT series can be reorganized, in what they named "renormalized PT" 
(RPT), in such a way as to obtain a better behaved expansion, valid in a larger range 
of scales. In particular, they showed that a certain quantity, namely the propagator 



defined in eq. (18), which measures the sensitivity of density and velocity perturbations 
to a variation in their initial conditions, can be computed exactly., i.e. at all orders 
in PT, in the large wavevector k (small scale) limit. CS finding prompted a certain 
number of independent investigations on possible ways to resum PT contributions at 
all orders, both for the propagator [221 ESI EH [25], and for the directly measurable 
PS[I11E61EIIEH1E91I3Q1I3II|32]. An alternative method, having the Zeldovich 
approximation as starting point, was proposed in [33l |3l]. The status of these methods 
to date can be summarized as follows: at 2; ~ 1 the PS can be computed at a few 
percent accuracy (in comparison with state of the art N-body simulations) in the BAO 
range of scales (0.05 ~ ~ 0.25h/Mpc), the accuracy degrading quite rapidly at higher 
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Figure 1. (left) The nonlinear PS from the N-body simulations of [35], divided by the 
no- wiggle PS of [36], plotted against the variable y defined in eq. ([l]). The color-code 
is the following: green for z = 0, black for z — 0.5, red for z — 1, and purple for z = 2. 
Each PS has been truncated at fc = 1 h/Mpc. (right) The relation between y and k at 
different redshifts. 



wave numbers (smaller scales) and smaller redshifts. Moreover, the computational time 
of these approaches, though much smaller than for N-body simulations, is still in the 
few hours range for a single PS, thereby making the implementations of MCMC's still 
quite problematic. 

In this paper we present a new computational scheme which overcomes the present 
limitations of resummed PT approaches in both respects: it greatly enlarges the range 
of scales in which it gives results accurate at the percent level, and it greatly reduces 
computational times. At z = 1 we can compute the PS up to A; ~ 1 h Mpc~^ at the 
percent level, in a time comparable to that of a simple 1-loop computation, namely, 
0(1) minutes. This opens the road both to parameter estimation via MCMCs, and to 
the extension of these methods from the BAO physics to weak lensing measurements. 

At the technical level, our main result is a resummation of the dominant PT 
corrections to the PS in the range of scales mentioned above. In the CS resummation 
for the propagator the effective expansion parameter in the large k limit turns out to 
be [21] 

y^e'ia^k, (1) 
where rj = log D{z)/D{zin) {D{z) being the linear growth factor), and 

with P^{q) the linear PS evaluated at the initial redshift, to be formally sent to infinity 
(in practice, in this paper we will use Zin = 100). Our starting point is the empirical 
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realization that, at wave numbers k larger than the BAO range, the nonlinear corrections 
to the PS take the form of a multiplicative function of the variable y defined in eq. ([!]). 
In fig. [l] we plot the ratio between the nonlinear PS from the N-body simulations of |35] 
at different redshifts, and the smooth linear PS by Eisenstein and Hu [36], as a function 
oi y. It is clear that, for y > 0(1), all the ratios are well approximated by an universal 
function of y. Therefore our goal is to identify, at each order in PT, the leading large 
y corrections which, once resummed at all orders, give the y-function plotted in fig. [T] 
In the range of scales considered in this paper {k ~ Ih Mpc-i) it is given, as a first 
approximation, by the simple analytic expression in eq. (56), which is obtained as a 
limit of the more refined eq. (53) discussed in sect. [6] 

We formulate a very efficient way to reorganize the PT expansion, based on 
evolution equations in time which are exact at all orders in PT. A similar approach 
was already presented in ref. [22] for the propagator, where it was used to reproduce the 
CS result, and also to include next-to-leading corrections. Here we derive the relevant 
equation for the PS and, after implementing an approximation analogous to the one 
leading to the CS result, we derive its behavior in the large k or, better, large y, limit. 
At the same time, at low k, the equation is solved by the 1-loop PS, thereby providing 
an interpolation between the two correct behaviors in the two extremes of the physically 
interesting range of scales. We should warn the reader here, as we will repeat again, that 
the range of scales where the eulerian approach is applicable is limited in the UV (small 
scales) by multi streaming, i.e. by small scale velocity dispersion, a physical effect which 
is absent from the approach even at the non-perturbative level. Nevertheless, the very 
large k limit or, more precisely, the limit in which all mode coupling is only with wave 
numbers q much lower than the scale k one is interested in, can be discussed exactly, and 
provides an useful starting point for the description of the intermediate scales, where 
the Eulerian approach is expected to work. 

The evolution equation contains two kinds of kernel functions. One is the same 
entering the equation for the propagator, and was already discussed in ref. [22]. A 
computation of this kernel function at 1-loop provides, by virtue of our evolution 
equation, the CS propagator at all loop order, as well as part of the leading corrections 
to the PS. The resummation of the leading corrections of the remaining part, which is 
responsible for mode-mode coupling, is the main result of this paper. 

As we have already anticipated, our combined resummation and interpolation 
procedure works remarkably well: its performance degrades only for (low) redshifts 
and (small) scales where the effect of multi-streaming, the intrinsic limit of eulerian 
PT, is known to become relevant [23 EB]- In other words, the comparison with high 
resolution N-body simulations presented in this work shows that, at the few percent 
level, our approach reproduces well the physics contained in the Euler-Poisson system, 
on which eulerian PT is founded. As a consequence, it provides the best starting base 
for methods aimed at going beyond the single-stream approximation, as that proposed 
in [2S|. 

The paper is organized as follows. In sec. [2] we review eulerian PT in the compact 
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form discussed in [211 |26], the diagrammatic approach to the PT expansion, and we 
recall the exact expressions for the fully nonlinear propagator and PS, which provide 
the starting points for their evolution equations. In sec. |3] we derive the evolution 
equations for the propagator and for the PS. In sec. |4] we show how, in the — )• 
limit, the PS equation is solved by the standard 1-loop expression for the PS, and in 
sec. [5] we discuss how a simple upgrading of the 1-loop expression provides an all-order 
resummation valid in the small k regime. Then, in sec. |6] we discuss the large k limit 
and discuss how to interpolate between this and the low k regime. In sec. [7] we present 
our numerical results and their comparison with N-body simulations. Finally, in sec. |8] 
we discuss our results and possible developments. 

In [Appendix A we discuss an alternative framework to do RPT, and in Appendix B 



we prove an useful relation valid in the large k regime. In Appendix C and Appendix D 



we give all the explicit formulas needed for the numerical implementation of the method. 
2. Nonlinear fluid equations and perturbation theory 

Eulerian PT ^ aims at solving the system of three fluid equations: continuity, Euler, 
Poisson. Starting from the case of Einstein-de Sitter cosmology, the equations can be 
written as follows 

^ + V-[(1 + Uv] = 0, 

OT 

^ +?/v + (v ■ V)v = -V0, 

OT 

V^(f)=^nUm, (3) 

where H = d log a /dr is the Hubble Parameter in conformal time, while (5m (x, r) and 
v(x, r) are the DM number- density fluctuation and the DM peculiar velocity field, 
respectively. 

Defining, as usual, the velocity divergence 6'(x, r) = V ■ v(x, r), and going to 
Fourier space, the equations in ^ can be expressed as 

95m(k,r) 



dr 



+ j rf^q(i^p5z)(k-q-p)a(q,p)6'(q,r)(5m(p,r) = 0, 



The nonlinearity and non-locality of the fluid equation are encoded in the two functions 



(p + q)q 0, ^ (p + q)^pq 



«q,p = 2 — , pq,p = — ^-^1 — , (5) 

which couple different modes of density and velocity fluctuations. 
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One can write Eqs. Q in a compact form [21]. First, we introduce the doublet 
(a = 1,2), given by 

MKv)\^ -n( S^iKv) \ (f.^ 

where the time variable has been replaced by the logarithm of the scale factor, 

1 " 

^ = log — , 

Qin being the scale factor at a conveniently remote epoch, in which all the relevant scales 
are well inside the linear regime. 

Then, we define a vertex function, 7abc(k, P,q) (a, c, = 1,2) whose only 
independent, non- vanishing, elements are 

7121 (k, p, q) = i 5^(k + p + q) a(p, q) , 

7222(k, p, q) = 5D(k + p + q)/3(p,q), (7) 

and 7121 (k, p, q) = 7112 (k, q, p). 

The two equations Q can now be rewritten in a compact form as 

dr, V2a(k, 1]) = -^lab ifbiK V) + e'^labciK -P, -q)Vb{P, V) V'c(q, ??), (§) 

where 

/ 1 -l\ 



(9) 



3 3 
V "2 2 / 

Repeated indices are summed over, and integration over momenta q and p is understood. 

To extend the validity of this approach to ACDM, we will reinterpret the variable 
rj as the logarithm of the linear growth factor of the growing mode, i.e. [U [T3], 

r] = \n{D/Din) , (10) 
and we redefine the field in Eq. ^ as 

with f = d\nD / d\n.a. As discussed in |14j, the above approximation is accurate at 
better than 1% level in the whole range of redshifts and scales we are interested in. 

If we consider the linear equations (obtained in the e^'jabc limit) we can define 
the linear retarded propagator as the operator giving the evolution of the field ipa from 

Vin to 1], 

'PaiK V) = 9ab{V, Vin)<fl{K Vin) , (12) 

where the "0" index stands for the linear approximation to the full solution. The linear 
propagator obeys the equation 

(Sabdri + ^ab)gbci'n, Vin) = ^ac^oiV ~ Vin)- (13) 
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propagator: -igabiv 



a power spectrum: P°ft(k, 77a, 77b) 



I 

A 



interaction vertex: -ie'^jabd^ay^b, kc 



Figure 2. The Feynman Rules for cosmological perturbation theory 
q k-q 



... k 




s k-q 




Figure 3. The two IPI functions Y^l^J and ^''^J at 1-loop, given expUcitly by eqs. (C.2) 



and (C.4) 



with causal boundary conditions. It is given explicitly by the following expression 



B + Ae-^/2{,-,') 



- ah 



div - v') 



with 9 the step-function, and 
1/32 



B 



5 \ 3 2 



and 



1 / 2-2 

-3 3 



(14) 



(15) 



The growing {ipa oc const.) and the decaying {{pa oc exp(— 5/2(7/ — rj'))) modes can be 
selected by considering initial fields ipa proportional to 



and Vr, 



1 

-3/2 



(16) 



respectively. 

In the following, we will be interested in the PS 

(V9a(k, T])ipb(k', 7]')) = ^^(k + k!)Pab{k] T], T]') 

and in the propagator, 

{J^^^^^)=^So{k + k)G.b{k■,v,v), 



(17) 
(18) 
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which gives the response of the field at time r/ to a change in the field at an earlier time 
rj' < rj. The linear approximations of the fully nonlinear quantities above are given by 
the linear PS 

Pab{k; V') = 9ac{V, Vin)9bd{V, Vtn)Pcdi.k] Vin, Vin) , (19) 



and by the linear propagator gabiViV') of ^q. (14), respectively. 

These quantities, as well as all higher order correlation functions, can be computed 
perturbatively by solving iteratively the nonlinear equations (|8|, supplemented by the 
initial conditions at the time rjin for the PS, the bispectrum, and, in principle, all higher 
correlation functions. In this paper, we will limit ourselves to the case of gaussian initial 
conditions, and therefore we will assume that the initial conditions are fully characterized 
by the PS only. 

A very convenient way to organize the perturbative expansion is to use the 
diagrammatic language discussed in [20] and in [26] (for details on the derivation of 
formalism used in this paper, the reader is referred to the latter). The building blocks 
are given in fig. [2| and are the linear propagator, the linear PS, and the interaction 
vertex of eq. ([T]). Propagators and PS can only be joined according to the rule encoded 
in the vertex, i.e., two continuous ends join with a dashed one. Notice, moreover, that 
according to eq. (19), a linear PS can be seen as the initial PS, represented by the empty 
square in fig. [2| sandwiched between two linear propagators. This will be useful in the 
following. 

The full nonlinear PS and propagator have the following structures [26] 



ds ds Gacik] T], s)Gbdik; r]\ s')^cdik; s, s') 



and 

Gab{k; V, V') = [9ba - ^ba] {k; 7], 7]') , 

where the last expression has to be interpreted in a formal sense, that is, 
Gab{k] r],r]') = gabiv - v') 

+ j ds ds' Qaciv - s)i:^d{k] s, s')gdb{s' - i) ^ 



(20) 
(21) 



(22) 



Notice that, both in (20) and in (21) and everywhere else unless explicitly stated 
otherwise, the time integrals run from — oo to +00. However, in practice, due to 
the causal properties of the retarded propagators (see eq. (14)), the range of time 
integrations turns out to be always finite. 

Eqs. (20) and (21) can be obtained either by analyzing the diagrammatic structure 
of the contributions at an arbitrarily high order, as in [20], or by functional methods, as 
in [2S]- They are exact relations {i.e. valid at all orders in PT), and are valid both for 
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gaussian and for non-gaussian initial conditions. Eq. (21) has been already exploited 



as the starting point of the derivation of the evolution equation for the propagator 
discussed in [22]. In the next section we will derive the evolution equation for the PS, 



starting from eq. (20). 



In the following, we will take the initial conditions for the PS, Pcd{k',riin,riin), to 
be given by the linear PS in the growing mode, P^{k)ucUd, which is equivalent to 
neglecting any nonlinear effect in the initial PS on the range of scale of interest. This 
is of course exact in the rjin — t- — oo limit. The linear PS in eq. (19) then becomes 



Pabi^jV^v') = UaUbP'^{k), where we have used the property of the linear propagator, 
9ac{'niV')'^c = which can be checked using eqs. (14) and (16 if] 



The crucial quantities entering eqs. (20 ) and (21 ) above are the 1-particle- irreducible 



(IPI) functions T^ab and where IPI means, as usual in Feynman diagrammatic 
language, that the diagrams contributing to these quantities cannot be separated into 
two disjoint parts by cutting a single PS or propagator line. T^abik] s, s') connects a 
dashed end at time s to a continuous one at time s' < s, whereas $ab(fc; s, s'), connects 
two continuous lines at s and s' with no time ordering. The lowest order contributions 
to these IPI functions are represented by the diagrams in fig. |3| and correspond to the 



explicit formulae given in Appendix C 



The information contained in the IPI functions T^ab and ^ab is completely equivalent 
to that encoded in the connected ones, namely the propagator and the PS. However, 
working with the IPI functions inside the time evolution equations to be introduced in 
next section presents some advantage, already exploited in |22]. Indeed, this approach 
provides a natural guide to interpolate between the — )■ and — )■ oo limits and, 
moreover, it is more suitable to be used for cosmology with time and/or scale dependent 
linear growth factors, which can be encoded in a time and scale dependence on the flab 
matrices ([9]), see [H]. 

The main focus of this paper will be on the nonlinear PS evaluated at equal times 
rj = rj' (hereafter, the PS), that is, 

Pab{k;v) = Pab{k;v,v)- (23) 



3. The evolution equations 

In order to derive the time-evolution of Pab{k', rj) we need that of Gab{k', rj, rj'), which can 
be obtained by using the property of the linear propagator 

drj gabiV, V') = Sab 5d(^? " V') " ^ac gcbiV, V') , (24) 

in eq. (g. We get (see P]), 

dr, Gab{k; T], 7]') = 5ab Sd{V - 7]') - Vtac Gcb{k] T], T]') + AGab{k; T], 7]') , (25) 

where 

AGab{k; r],ri') = j ds T.ad{k; r],s')Gdb{k; s',r]'). (26) 
I Remember that our fluctuations, defined in eq. ml, are constant at tlie linear level. 
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Notice that the terms containing the Dirac deha in (24) and (25) where not written in 



[22] . because in that paper we always considered propagators at different times. Here, 
on the other hand, those terms are important when the propagator is inside a time 



integral, as in the second hne of eq. (20). 



Taking the rj derivative of Pab{k', rj) as given by (20), we get 
dr^ Pab{k; rj) = -Qac Pcb{k; r]) - Qbc Pac{k; r]) 
+ ( AGac(fc; ?7, Tlin)Ghd{k] ?7, r]in) + Gac{k] r], riin)AGM{k; r], T]in) ) P^{k)ucUd 



+ / ds' [^ac{k;ri,s')Gbc{k;ri,s') + Gac{k; r^, s')$cfe(/^; s', r/)] 



+ j dsds ^cdik;s,s')(^AGac{k;'n,s)Gbd{k; rj^s 
+ Gacik;T],s)AGMik; r],s' 



which is of course still an exact equation. 

In [22], we showed that in the large k limit AGac takes the form 

AGac{k; T], s) ~ Ha{k; i], s) G^dk; r], s) , 

where 



(27) 



(28) 



with S^^-* the 1-loop approximation to the full T^ad, see eqs. (Kill) and (|C.2[), and the 



boldface index indicates that we are not summing over that index even if it is repeated. 



The same factorization, eq. (28), holds at A; — > where the 1-loop PS is obtained from 



the solution of the evolution equation. Therefore, in [22|, a natural way to interpolate 



between the two extreme limits has been identified in using eq. (28) also for intermediate 
fc's. 



We will denote the solution of eq. (25) with the approximation (28) as Gab{k; rj, rj') 



which we will use in (27). As it was discussed thoroughly in [22], the solution Gab{k; rj, rj') 



is exact both in the low and in the large k limits. At low k it reproduces the 1-loop 
propagator 

Gabik; r],r]') gabiv,v') 
whereas at large k it reproduces the exact result obtained by CS in [21] 



(29) 



Gab{k; T], T]') Gil (k; r], rj') = Qabij] - v') exp 



2 2 ' 



(30) 



(for k — 7- oo) , 

where is defined in eq. ([2]), and "eifc" stands for the eikonal limit, in which the above 
expression is exact [ID], SI], see Sect. [6) 
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4. Recovering 1-loop 



As a first attaclc to a practical and fast solution of eq. (27), we identify the limit in 



which it reproduces the 1-loop result. It is obtained by using eq. (28) in (27) and by 
setting 



Gab gab , ^ab $11^ , (31) 

in the second and third lines, where is the 1-loop approximation to ^aft- Moreover, 
the fourth and fifth lines, containing at least 2-loop order quantities, are consistently 
neglected. Thus, we get 

d,P^c^{k-r,) = -QacP^\k;v)-nbcP!.l\k;v) 

+ P°(A;) (^Ha{k; r], T]in)ub + Hb{k; r], r]in)u^ 

ds (^^l^Jik;ri,s)gbdiv,s) + gadiv,s)^M\k'^s,ri)^ . 

(32) 

The solution of the equation above exactly reproduces the 1-loop PS, 



UaUb 



+ j ds{gae{ri,s)ub + gbe{il,s)ua) He{k;s,riin) 
+ j ds ds'gac{v,s)gbd{v, s')^^^{k;s,s') . 



(33) 



as can be checked directly by taking the ?7-derivative of the expression above. 



5. The small k limit 

The approximations leading to the 1-loop result suggest the first step to take to 
go beyond, and to obtain a first resummation, containing infinite orders of the PT 



expansion. Indeed, by promoting the linear PS appearing in eq. (32) to the nonlinear 
and time-depedent one, we get 

d, PS>\k; - Pi\k; - Qbc P,f{k; r^) 

+ H^{k- r],r],^)P^\k;r]) + Hy,{k- r],r]i^) P^^{k-r]) 

+ I ds i^'ahk- T], s)gbd{v, s) + gadiv, s)'^'!!^^; s, v)] ■ 



(34) 

To understand what this approximation corresponds to, it is instructive to set to zero 
the third line of the above equation. Then, the equation can be solved exactly, to get 

^a1^(^;^)|$(i)=o = Gacik; r],r]in)Gbd{k] r],r]in)ucUdP°{k) . (35) 
The above expression contains infinite orders in PT, and corresponds to the linear 
PS multiplied by two renormalized propagators. By turning ^^^^j on we are adding 
perturbatively the effect of mode-mode coupling. 
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k 




q 



A B C 

Figure 4. The 2-loop contributions to ^ab giving leading contributions in the large k 
limit, with no vertex renormalization. 




q 



Figure 5. The dominant lowest order contribution to $afc with a vertex 
renormalization. 



A step further, which does not increase the computing time too much, consists in 
improving the approximation in eq. (31), by using the 1-loop approximation for Gab, 



instead of the hnear one, in the second hne of eq. ([34j), to get 

\k; r^)= - nacPi\k;r^) - n^cP^fik; r/) 



dr,P, 



R2f 
ab 



tdR2 
ab 



+ I ds[<l>^^J{k;r^,s)Gl'>{k;v,s) + G^:>{k-r^,s)<i>i;,\k-s,v)] 



(1)^ 



(36) 



where the 1-loop propagator has been given in eq. (29). 

Numerical results for the 1-loop approximation of eq. (33) and for the small k 
resummation discussed above, P^^, will be presented in Sect. 7j 
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m,pairings ^ 



# : renormalized vertex 




+ 




Figure 6. The diagrams contributing to ^"^'^{k\ri^s')G^^^{k\ri,s'). Thick lines and 
dots represent, respectively, renormalized propagators and vertices in the eikonal limit. 
The sum over m and over all possible pairings of soft lines can be expressed as the two 



diagrams at the RHS. The black PS is given by eq. (46) 



6. The large k limit 



6.1. The eikonal limit 



We will now derive tlie large-fc approximation of tlie equation for tlie PS, extending 
the derivation of ref. |22] for the function T,ab{k;ri,ri') to the mode-coupling function, 
^ac{k',f],'ri'). In ref. |22], the inclusion of all the dominant configurations in the large k 



limit allowed us to write eq. (26) as (28). In the present case, we will compute taking 



into account the diagram on the right of fig. [3} in the k ^ q limit, and all those diagrams 
at all orders, that can be obtained from that one by attaching "soft" PS's, with momenta 
q k in all possible ways to the "hard" lines carrying momenta of 0{k), both inside the 
loop and as corrections to the vertices. For instance, at 2-loop order, the diagrams we 
will take into account are those of figs. |4] and |5} We will denote the resulting ^ac{k; rj, rj') 
as ^^^{k; rj, rj'), since it is computed in the "eikonal" limit in which the wavenumber we 
are interested in, k, is much larger than all the wave numbers it is effectively coupled 



to, see refs. [lOl HI]. In Appendix A we discuss a framework to do renormalized FT, 
that we call eRFT, defined as a loop expansion around this extreme situation, where 
loop corrections progressively take into account the effects from intermediate modes, 
and restore the full momentum dependence at small and intermediate scales. 



As we will show in [Appendix B[ in the eikonal limit the n— th order contributions 



(in standard FT) to the last three lines of eq. (27) can be written as 



n-l 



dsds' J2{^GZ''^'\k-,v,sW:,\k-,s,s')Glf{k-, v,s')] 
1=1 



(n-l) 
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feik 



eik I 



\eik I 



in) 



n-1 



dsds' Y.{Hfik; V.mn)Gt^'^'''\k- r^,s)[^-'{k;s,s')Glf{k; v,s')] 
1=1 



/\]in-l) 



[GZ\k; V,s)^:,\k;s,s')f-^^Gf/^-'\k-, v,s')H-\k; v,Vr, 



+ 



ds' 



<l>Z''\k- V, s')GZ\k- r^, s') + GZ\k- r^, s')<!>f'''{k- s', r,) 



(37) 



where Hf^{k] r], r]in) = J^^^ ds'Yfae'^^\k] r], s')ue = -Ua k'^al e'^{e'^ - e'''"). 

The first two hnes at the RHS of (37), summed over n and combined with the 
second hne of (27), give 



yeik 



eik I 



(3J 



where we have used the fact that P^l'^ has the structure given in eq. (20), with the full 
Gab and replaced by their eikonal limits. 

The functions ^^b '^ ^'^ab '^ ^'^^ obtained from the full by taking into 
account only those contributions having a tree level, i.e. not-renormalized, vertex at 
the end corresponding to the index "a" and "b" , respectively. The first term in the last 
line of eq. (37) can be schematically represented as in the LHS of fig. |6} There, the 
horizontal thick lines represent eikonal propagators as defined in eq. (30), carrying the 
hard momentum fc, whereas the thin lines represent soft propagators, which can be well 
approximated by the linear gab- 

The m lines on the left have to be joined to the m ones on the right in all possible 
ways, by means of m soft PS's. Notice that the number m in the figure does not 
correspond to the order in standard PT, since the thick lines and the thick dots already 
include infinite orders in standard PT. In particular, the thick dots represent fully 
renormalized vertices in the eikonal limit. These vertices are obtained by correcting the 
tree level expression 



1 k ■ q 

z q 



(for k > q) , (39) 



by including all possible soft PS insertion on the hard propagator lines. At 1-loop, the 
expression above is corrected by the three diagrams in fig. [7| which, again in the eikonal 
limit, sum up to. 



2 2 ' 



k'a^ 



Ik- q 



Q{ri — s)6(s — rj') , 



(40) 



where we recognize, at the first factor, the 1-loop contribution to G^^{k;ri,ri'). 
Considering all loop orders, and summing up, one finds 

ds, ds2 G'J^{k; T], si)r^f^(k, -q, q - k; Si, s, S2)u, G^f (Iq - k|; Sa, r/') 
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n' 



._A.J.._' 




Figure 7. The one- loop correction to eq. ( 39 ) 



s)Q{s - i) (for k > q) , (41) 



where r^*^g is the fully renormalized vertex in the eikonal limit, represented as a thick 
dot in fig. [gJ Notice that the leftmost (rightmost, if is considered) vertex is not 
thick, i.e. renormalized, but it is given by the tree level 7af,c of eq. ([T]). 

Therefore, in the eikonal limit, the time integrations on the insertion points of the 
soft legs via a full vertex, i.e. those corresponding to the thick dots in fig. [6| factorize 
from the time integrations of the loops correcting QabijIiV') to C^abi.^jVyV')y can be 
performed independently, provided the causal ordering enforced by the theta-functions 



in eq. (41) is respected. 



The contributions obtained by taking into account the ml possible pairings between 
the soft lines by the m soft PS's in fig. |6} and then by integrating the m — 1 full vertices 
on the left and the m ones on the right from r]in to t], give 



ml [m — 1)! 

which, summed on m from 1 to oo gives the simple result 

ds [K!''ik;V,s)Glf{k; V,s) + GZ\k; ry, .)<'^(A;; r/)] 



(42) 



= 2P'^{k)uaUbk^ale 
where we have taken the rjin — ?■ —oo limit. 



(43) 



In the eikonal limit, the evolution equation (27) can then be rewritten as 

dr,P:t{k;v) = -nacP:tik;v)-^bcP:^{k;v) 



+ Hf{k- r^, r^,„) (fc; r^) + H^^'{k; v, mn) Pl\k- r^) 



yeik 



eik I 



+ 



ds [^f,\k-rt,s)Gt,'{k- v,s) + GZ'{k- v,s)<!>-t"{k-s,v)] 



yeik 



eik I 



ik.Ri 



O peifc 
,2 2 



{k-r^)~n,,P:f{k-r^) 



-2k:^ate^^P:i\t,r^) 
+ 2P\k)uaUbk^ay\ 



(44) 



where we have consistently taken Ha — ?■ H'^ also at the second line of (27). Notice the 



formal similarity between eq. (44) and the equations we derived for the low-/c regime. 



eqs. (34), (36), suggesting that, as for the propagator [22], an equation of the same 



structure could be used to interpolate between the two extremal regimes. 
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The equation above has the attractor solution 



->eik I 



.,,\k-ri) = P\k)u,u,, (45) 

as can be verified directly, by using fiacre = (see eqs. ([9]) and (16)). This result can 
be generalized to the two-point correlator at non-equal times ?7 and rj . The relevant 
diagrams to compute this quantity are obtained by attaching an eikonal propagator 
(connecting the time 77 to r/') to the left end of the LHS of fig. [6] and by replacing the 
leftmost tree level vertex by a thick one. Proceeding analogously to what we did above, 
one gets 



P:i^{k-r^,r^') = P%k)uaU,exp 



22 ' 



-k'a. 



(46) 



which reduces to (45) for rj = rj'. 

The equality between the fully nonlinear PS, in the eikonal limit, and the linear 
on^ can be understood physically as due to the fact that all the soft modes, in the 
extreme q <^ k limit, are seen, at the short scale set by /c, as zero modes. Since, by 
Galileian invariance, equal-time correlators should be independent on zero modes [55] , 
all the nonlinear corrections cancel out in this limit (see also |41j). 

In order to outline the interpolation procedure between the large and the small k 
regimes, it is useful to notice that the contributions leading to eq. (43), once summed 
over m, can also be represented by the two diagrams on the RHS of fig. |6| where the 
black PS now is the fully renormalized two-point correlator in the eikonal limit, eq. (46), 
and all the vertices besides the leftmost one are fully renormalized (in the eikonal limit). 
The first diagram comes from all pairings in which the vertex in rj is paired to that in 
s' the second one from all the other ones. 



6.2. The intermediate k regime 

The eikonal limit provides a conceptually interesting benchmark, but it does not 
adequately describe the large k limit. One one hand, in practical situations such as 
ACDM, the separation of scales is not so abrupt as to motivate an effectively g — )■ 
limit for the soft modes when the hard one is k = 0(1) hMpc"^ On the other 
hand, considering arbitrarily large momenta is meaningless because of multi streaming 
effects, which are totally neglected in the PT framework, but are certainly relevant for 
A; ~ 1 h Mpc"^ and higher [371 [38l [39] . 

In order to recover the proper momentum dependence for the nonlinear propagator 



at lower k, in [22] we used the 1-loop expression for Ha in the evolution equation (25). 



The equation so obtained provides an interpolation between the large k limit, in which 

-<eik 



^eik generates the all-order resummed propagator G^^bi ^^"^ 1-loop standard PT at 



small k. 

§ We are considering a single fluid with initial conditions in the adiabatic growing mode, for the case 
of isocurvature modes, see [101 HT] , 
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: renormalized vertex 



Figure 8. The diagram for E™^'^(fc; ry, s')Gj;f (fc; c', ?]') 



In this paper we will follow an analogous path to get an interpolation procedure 
for the nonlinear PS. First, we will assume that the evolution equation has the same 



and , appearing in those equations. As for Ha we will use 



structure as eqs. (34), (36), and (44), for any values of fc, and use appropriate values for 



the functions Ha, Gab 

the complete 1-loop expression, as we did in [22]. In principle Gab, appearing at the third 



line of the evolution equation, should be taken equal to Gab introduced in eqs. (29), (30) 



however we checked numerically that using the more numerically convenient G'^ab gives 
undistinguishable results. Then, it remains to be discussed the interpolating expression 



for the mode coupling function between the low-k value, $1^^ and the eikonal 

limit. Before doing that, as a warm-up, we consider again the evolution equation for 



(1) 



the propagator, eq. (25). 



In the evolution equation for the propagator, a time-integrated product between 
and Gcbi see eq. (26), appears, which in the eikonal limit takes the factorized expression 
in eq. (28). This quantity can be represented by just the 1-loop diagram of fig. Isl 



where the hard propagators and the rightmost vertex are the fully renormalized, eikonal 
ones. Notice that, being the vertex at s' a renormalized one, the two hard propagators 



combine in a single one, and the integration over s' factorizes (see eq. (41)), leading to 
eq. (28) with Ha = H^^^. Now, to get the correct — )■ limit from this expression one 



has to replace the fully renormalized vertex with a tree level one, obtaining in this case 
eq. (28) plus subleading 0{k^a^) corrections. This is the same result that one would 



obtain in eRPT at 1-loop, see Appendix A Therefore, the interpolation between large 
and small k can be entirely loaded on the shoulders of vertex renormalization. Indeed, 
at the scales of interest, vertex renormalization is not fully at work as to give the eikonal 



expression in eq. (41). This is shown in fig. ^ where we plot eq. (41 ) (orange line), the 



corresponding quantity in which the thick vertex has been replaced by a tree-level vertex 
(red line), 

r^, s)e^7-ice(k, -q, q - k)w, G^f (A;; s, V) , (47) 
and the improved interaction term 

Ikq 



eq. (41) 



Geik 
ad 



{k;ri 



Geik 
db 



{k;s,r]') 
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yq=0.5, x=0.1 



eq. (41) 

■ eq. (47) 

■ eq. (48) 




yq=0.5, x=0.9 



eq. (41) 
eq. (47) 
eq. (48) 




Figure 9. The improved vertex interaction of eq. (48) (purple line), the full eikonal 



expression of eq. ( 41 ) (orange line) , and the tree-level vertex interaction of eq. ( 47 ) (red 
line), for different values of and x = k ■ q^/kq. Notice that, at z — 0, i/q = 0.1, 0.5, 1 
correspond to q = 0.016,0.08,0.16 hMpc^^, respectively. 



+ GZ{k- T], s)e'jdce(k, -q, q - k)^^ GX{k- s, r,') 



(4J 



(purple line) obtained by replacing the tree-level vertex contribution to (41) in the 
eikonal limit with the tree level vertex in its complete form 7(ice(k, — q, q — k). Of course 



eq. (48) goes to eq. (41 ) in the eikonal limit, but the question here is to understand how 



far from this limit it is at the scales of interest in the paper. 

All the three expressions are computed setting the index a = 1, and have been 
contracted with ui, and integrated in s from rj' — )■ — oo to rj. Moreover, we have divided 
by exp(r7). The resulting quantities depend on yk = ka^e^, yg = qayC^, and x = k-q/kq. 
The plots in fig. |9] clearly show that the contribution from vertex renormalization (in the 
eikonal limit) is subdominant for yk up to 6, corresponding, at 2; = to 1 hMpc~^. Only 
for very small values of the soft momentum, the eikonal limit for the vertex interaction 
is approached. 

We are therefore motivated to extend the procedure already discussed above from 

f ( R) 

the to the , namely, neglecting vertex renormalization. This is further justified 
by the fact that the soft momentum contributing to is typically larger than those 
dominating the T,ab, due to the presence of a P^{q) in the loop. The prescription of 
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using tree level vertices is also motivated by eRPT, as discussed in Appendix A 

Looking at the RHS of fig. [6| we have two contributions. If we replace the full 
vertex in the first diagram with a tree-level one we get the 1-loop expression in eRPT 



for of eq. (A. 20), which in the k ^ q limit can be approximated as 

\2" 



e_RPr,{l) 
ab 



{k;s,s') ~ <l>il^(A;;s,s')exp 



-k'a' 



e — e 



(49) 



where ^^^jj {k; s, s') is computed in 1-loop standard FT. As for the second contribution 
we also replace the rightmost vertex with a tree-level one. In order to avoid 2-loop 
integrations, we further simplify this contribution by taking the tree level vertices in 
their limiting form, eq. (39), and take P°(|k — q|) ~ P^{k) inside the loops. We get 

\2" 



(feV^e'^+^'j P\k)uaUbexp 



(50) 



Since both contributions, eqs. (49) and (50), have tree level vertices at both ends, there 



is no difference between <l>^^ and at this level. The latter will therefore both be 



approximated by the sum of eqs. (49) and (50), namely 



^abik;s,s') = e 2 



+ ( k'a' e 



P\k)UaU, 



.(51) 



At small /c, the above expression explicitly goes to + 0(A;^ 



The second term 

inside parentheses comes from the second diagram at the RHS of fig. ^ and it is 
crucial in order to recover the eikonal limit at large fc, once vertex renormalization for 
the rightmost vertex is restored, as in fig. [6j However, it should be switched-off for small 
momentum values, because in that range it becomes degenerate with 2-loop diagrams 
not included in this resummation, which would give other 0{k'^) contributions. In order 
to avoid taking into account these new contributions, we suppress the second term in 



eq. (|51|) at small k multiplying it by a filter function of the form 

F{k) 



1 + {k/k) 



T\4 ' 



(52) 



see also eq. (D.7), where the power 4 ensures a rapid enough switch-off of the filter at 
low k. 



As we will see in next section, see fig. 12, the filter becomes irrelevant at large k 



and for redshifts ~ 1. Therefore we fix /c at z = (where the filter effect is maximal) 



by taking it equal to the scale k at which the two terms of eq. (51 ) are equal. Choosing 



a smaller value would suppress the 1-loop contribution to eq. (51) which is shown to 



successfully reproduce the BAO scales by our discussion of the small k range (see approx. 
"R2" introduced in sect. [5]). On the other hand, choosing a larger A;, would delay too 
much (in k) the effect of the large k contribution. This procedure gives the value 
k = 0.2 hMpc~^, that we will use explicitly in our numerical computation. In the next 



section we will also show results in which no filter function is used in eq. (51 ), see fig. 12 



and in Appendix D we will show results obtained by varying the value of the filter w.r.t. 
the optimal value. 
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In summary, in the following we will solve the evolution equation 
dr, Pab{k; ri)= - Vtac Pcb{k; 77) - ^bc Pac{k; 1]) 

+ H^{k; 7], Tjin) Pe,b{k; v) + Hi,ik; r], rjin) Pah{k; v) 

+ I ds [^ad{k; r/, s)Glf{k; r/, s) + Gl'^{k; r/, s)^db{k; s, r])] 



(53) 

where Ha{k; rj^rjin) is given in eq. (28) and ^ad{k]rj,s) is given in eq.(51) with the filter 



function (|52|) inserted as described above. Explicit expressions for these functions are 
given in 



Appendix D 



We reiterate that the above equation interpolates between the small k and the large 
k limits discussed in the previous sections. For small /c, we get the approximation R2 
discussed in Sect. ^[[ 

The large fc, eikonal limit, is recovered by replacing back the full vertex r^*^,^ for the 
tree one '-)abc at the rightmost (the leftmost one, if is considered) vertices in fig. ^ 



This operation turns the last line of eq. (53) into the eikonal limit of eq. (44). However 



since, as we discussed above, the eikonal limit is not realized in practice, even at very 
large k, we will not introduce an explicitly interpolating procedure between these two 



expressions, but will use eq. (53) for all values of /c, using the comparison with N-body 



simulations to assess the range of validity of the approximation. An alternative, more 
systematic but also more computing time consuming approach, would be to consider 



the improved vertex in (48) and/or higher orders in eRPT. We leave the investigation 
of this way for future work. 



It is therefore instructive to study the very large k limit of eq. (53). Using the large 



k limit for the first term of eq. (51) 



KT^^^'\k;s.s) ^ [k^ale'^^') P{k)uaU, exp 



-k^al 



(54) 



the last line of eq. (53) can be integrated analytically in s, to give 

UaU,P\k) [y\e-y' -l) + ^y{l+y^)Y.ii{y)l where ?/ = eUa, , (55) 

and "Erf is the error function, Erf(?/) = 27r^^/^ dte~^^. At large y, the above 
expression goes as P^{k)y/7iy^{l — l/yy/n + 1/y^)- Considering that Ha and Hf,, in the 
same regime, go as — y^, we realize that in the large-y limit the differential equation in 



(|53|) has an attractor solution, given by the simple formula 

1 1 

+ ^ 



TC 



Pab{k; rj) ^y 



1 - 



UaUbP\k) 



for large y , (56) 



yv^ y 

which exhibits the pure y dependence for the nonlinear to linear PS ratio we found for 
the N-body simulations, see fig. [Tj 

In sect, [rj we will present results for the solutions of eq. (53) as well as for the 



analytic approximation (56). 



II Actually, to get exactly R2 in the small k limit, we should use G instead of G"^"^ at the last line 
of eq. (53 1, however the two choices give sub percent differences, while the use of G^^^ is much more 



convenient numerically. 
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Figure 10. PS's normalized to tlie no- wiggle PS of [35] at different redshifts. The 
color code is the following: the grey dash-dotted line is linear FT, the green dotted line 
is 1-loop PT, the blue dashed line corresponds to the solution of the evolution equation 
in the small k approximation "R2" of eq. ( 36 ) , the red solid line is the solution of the 
evolution equation ( 53 1 , and dots with error-bars are the N-body results of Sato and 
Matsubara [3S]. Also shown (thin grey line) is the output of the Coyote interpolator 
of refs. [HilliS]. 



7. Numerical results: comparison with N-body simulations 

The aim of this section is to test numerically the goodness of our approach. We 
compared the nonlinear PS predictions against the data coming from high accuracy N- 



Nonlinear Power Spectrum from Resummed Perturbation Theory: a Leap Beyond the BAO Scale22 





Figure 11. The relative difference between various approximations discussed in the 
text and the N-body simulations of Sato and Matsubara [3S]. Same color-code as in 
fig.[TOj The thin grey line is the comparison with the output of the Coyote interpolator. 

body simulations designed to predict the nonlinear PS at the percent accuracy around 
the BAO range of scales. The initial linear PS for the simulations we considered was 
obtained from the CAMB public code [H] and any primordial non-Gaussianity was 
neglected. Accordingly, we solved our evolution equations taking the initial PS from 
CAMB at Zin = 100, where the gravitational clustering is fully linear on the scales of 
interest. 

We provide plots of the comparison between our approach and the N-body 
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simulations presented by Sato and Matsubara in [35]. They considered a ACDM 
cosmology with the following parameters: Qm = 0.265, Qi^h"^ = 0.0226, h = 0.71, 
n = 0.963 and as = 0.8. In order to check the possible dependence of our results on 
the cosmology or on the N-body simulations, we also compared our results with the 
independent set of simulations produced by Carlson, White and Padmanabhan in 
for a different ACDM cosmology, and with the results from the cosmic Coyote emulator 

Wli US] . The latter is an interpolator built from a suite of many different N-body 
simulations which provides nonlinear PS's for different cosmological parameters chosen 
inside a certain range and for redshifts z < 1. In all these cases our comparison tests 
worked at the same quality level. 

In fig. 10 we plot the nonlinear PS computed in different approximations, divided 
by the smooth linear PS given in [SB]. The blue dashed line corresponds to the small 
k approximation of eq. (36), while the red sohd line is obtained by using eq. (53). The 
results obtained in linear PT (grey dash-dotted) and 1-loop PT (been dotted) are also 
shown. Relative differences with respect to the N-body results are given in fig. 11 with 
the same color-code. 

We notice that our evolution equation in the small k limit of eq. (36), i.e. with 
no-resummation of the contributions to ^ab beyond 1-loop order, is able to reproduce 
the nonlinear PS at the percent level in the BAO range. There are a couple of outlier 
points in the linear region k < 0.1 h/Mpc which can be ascribed to a fluctuation in 
the N-body simulations of [35] from a comparison with the Coyote emulator (thin grey 
line in fig. 10, 11) and with the 1-loop approximations, which works well in this region. 
Concerning the resummed and interpolated solution, i.e. the red line, we see that it 
performs at the 1% level in the BAO region at any redshift, including z = 0, where, due 
to the larger amount of nonlinearity, the last peak is at about 0.15 h/Mpc. At redshifts 
z > 0.5 it performs as the Coyote interpolator up to A; ~ 0.35 h/Mpc, namely, still in 
agreement with the N-body results of [35j at the 1% level. 

In fig. [121 we plot, in an extended k range, the results obtained as in the red line 



of fig. 10 with (thick red) and without (thin red) the filter function of eq. (52). We see 
that, as expected, the effect of the filter vanishes at large k (where the filter goes to 
unity) and at increasing redshift, where the relative weight of the second term in eq. (51 ) 
compared to the first one becomes less relevant. We also plot (thick black dashed line) 
the result obtained by using the large k expression of eq. (54) in eq. (51). We also show 
the analytic expression for the attractor solution found in eq. (56) (thin black dashed 
line) . 

The comparison with the nonlinear PS from the N-body simulations at large k 
is given in fig. 13 In order to gauge the performance of the N-body results at 



large k, Sato and Matsubara performed runs with two different volumes, a large one 
{Lfjox = lOOOMpc/h) and a small one {Lbox = 500Mpc/h), which we plot with black 
and blue points, respectively. We notice that the two sets of N-body data are practically 
overlapped for z < 1, where we also checked that they agree with the Coyote emulator, 
but diverge significantly at z > 2 for k ~ 0.8 h/Mpc. These transient effects therefore 
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Figure 12. The PS at large k from the evolution equation in various approximations 
for the mode-mode coupling part in the third line of eq. (53) and at different redshifts. 
The thick black dashed line corresponds to using eq. (55), the thin red solid line to 
the improved ^ab of eq. (51) with no filter imposed in the second term, and the thick 
red solid line to using eq. (51 1 with the power-law filter function of eq. (52 1 with 
k = 0.2h/Mpc (see Appendix D for explicit formulae). Also shown are the large y 
attractor solution of eq. ( 56 ) (thin black dashed line) , and linear PT (grey dash-dotted 
line) . 
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Figure 13. Comparison between the evolution equation discussed in this paper (red 
sohd line, see Appendix D for explicit formulae), and the N-body simulations of Sato 
and Matsubara black dots are for the large volume simulation {L — lOOOMpc/h), 
blue ones for the small volume one (L = 500Mpc/h). Also shown are linear PT (grey 
dash-dotted line) and 1-loop PT (green dotted line). 
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prevent us from considering the comparison for higher redshifts and scales, however, the 
trend from z = up to 2; = 2 clearly shows a progressive improvement, as it should. 
Quantitatively, at 2; = 1 we measure an agreement at the 1% level between our results 
and the N-body simulations up to A; = 0.8 h/Mpc and at 2% on the same range of scales 
for z = 0.5. 

8. Discussion and conclusions 

The scheme discussed in this paper presents three main advantages with respect to 
alternative approaches to the nonlinear PS in the ~ 1 hMpc~^ scale, including N-body 
simulations: the accuracy (already discussed in the previous section), the computational 
speed, and the range of cosmologies that can be dealt with. The computational time is 
the same as that required by a standard PT computation at 1-loop. Indeed, concerning 
momentum integrations, to obtain the three independent components of the PS, one 



needs to perform two one-dimensional integrals, eqs. (D.l ) and (D.2), and the three two- 



dimensional ones of eq. (C.4). Exactly the same integrals enter the 1-loop expression for 
the PS, eq. (33). In addition, the solution of our evolution equation requires only one 
more integration in rj, which however takes a time of the order of a second. Therefore, 
any point in k takes at most a few seconds to be evaluated. 

The extension from ACDM cosmology to more general ones is also greatly eased in 
this approach. In principle, all the cosmologies in which the fluid equations of eq. ^ 
are modified only in the linear terms, can be taken into account by a modification of 
the ilab matrices in (|9]), as it was discussed in [H]. Cosmologies of this type include, 
for instance, those with massive neutrinos |l6l STJ HH] , modifications of gravity of the 
scalar-tensor/f(R) type |49], or Dark Energy models with a non-relativistic sound speed 
[18] . Non-gaussian initial conditions can also be taken into account by the inclusion of 
new vertices in the diagrammatic rules of fig. |2] and their impact on the propagator and 
on the PS can be analyzed [50l [51] . 

Evolution equations for the nonlinear PS have been proposed also before [521 
[27[ dU [S3]- The main step forward provided by the present analysis is, besides the 
computational speed, the improved treatment of the large k limit, in particular for the 
mode-mode coupling term. Indeed, while the gaussian damping of the contribution 



proportional to the linear PS, namely the first term in eq. (20), is reproduced 
satisfactorily well by most resummed PT methods, the remaining part, that containing 
the mode-mixing function is usually included only at lowest order in PT, adding at 
most a few corrections. As a result, the improvement of these methods over the 1-loop 
results was not able to go beyond the BAO scale, as it was discussed in detail for the 
Time Renormalization Group approach of ref . [H] in [Ml 132] . 

A different strategy to resum PT corrections at all orders is provided by the 
approach of [21]. In this context, it was shown that the exact CS resummation for 
the propagator holds also for the larger class of multipoint propagators. The nonlinear 
PS can be expressed as an infinite sum over squares of multi-point propagators, however 
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the kinematical regime for which the exact resummation was proved is different from 
the one considered here, namely, the one in which only one internal PS is at large k, the 
other ones being soft. Our results, and their good comparison with N-body simulations, 
indicate that the relevant kinematical configurations in the intermediate k regime are 
still given by those loops in which most of the external momentum is carried by a single 
internal leg, and that the dominant diagrams are still of the chain type, with some 



amount of vertex renormalization neglected, as detailed in subsection 6.2 



Our results depend on the fact that ACDM-like PS's are well behaved as — )■ 
and have an intrinsic length scale which is of the order of a^e"^. This makes a separation 
between 'hard' and 'soft' modes possible, and, due to the k/q enhancement of the vertex 
functions, singles out the class of chain diagrams as the dominant one at large k. A very 
different approach would be needed to take into account, for instance, pure power- law 
PS's. In this case the dominant mode coupling at a scale k would involve modes of 0{k) 
as well. Fortunately such cases are nowadays only of academic interest. 

The approach presented in this paper suffers from an intrinsic physical limitation, 
namely, the neglect of velocity dispersion and all higher order moments of the particle 
distribution function, which is at the basis of the derivation of eqs. ^ from the Vlasov 
equation. This "single stream approximation" is known to hold at large scales and high 
redshifts, but it was estimated to fail at the percent level in the BAO range at 2; — )■ 
[371 EH]. The comparison of Sect. [T] between our results and N-body simulations exhibits 
the same trend. It will be interesting to investigate further the origin of this excellent 
agreement. 

The access to the multi-stream regime is precluded by construction to all forms of 
(improved) PT, including this one. A way to incorporate such effects in semi- analytical 
methods was recently proposed in [39] , in which the feeding of the multi stream at small 
scales on the more perturbative intermediate scales was described in terms of effective 
source terms. The inclusion of such effects in the present approach will be studied 
elsewhere. 

For the time being, the increase of the maximum k at which the nonlinear PS can be 
computed reliably, provided by our approach, opens the way to interesting cosmological 
applications, allowing a tighter extraction of cosmological parameters from the LSS 
relevant for topics such as the measurement of the acoustic scale from BAO's, the limit 
on the neutrino mass scale, and possibly, cosmic shear. 
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Appendix A. 

In this appendix we will discuss a scheme to do renormalized cosmological perturbation 
theory in a way that closely resembles the RPT of [201 HI] but gives a systematic way 
to compute higher loop orders. In particular, we will stress the crucial role of counter 
terms, which should be added to the diagrammatic rules in order to avoid overcounting 
of the contributions. Since the 'tree level' in this framework coincides, in the large 
momentum limit, with the eikonal approximation of [lOlHT] we will refer to this scheme 
es eRPT. 

The root of RPT as defined in [201 EI] is to use, as tree propagator, the resummed 
one i.e. the one decaying as exp [ — k'^a'^{e^ — e'''')^/2] at large k, instead of the linear 
one, Qabiv^v')- Since the low-fc limit of the resummed propagator does not match with 
1-loop standard PT, an interpolation procedure is usually advocated in order to have a 
well behaved tree propagator, going to the 1-loop (in standard PD) one as A; — )■ 0. 

An alternative approach is to define as the tree level propagator, for all values of 
k, the resummed propagator 

G2\k;v,v') = 9ativ,v') exp [ - k' a'^-^—f^] , (A.l) 
and as tree level PS, 

P^tik; V, V') = P'{k)uau, exp [ - k^al ^""' ~f ] , (A.2) 

where P^{k) is the linear PS in standard PT and we have assumed growing mode (and 
adiabatic) initial conditions. The above quantities are obtained starting from the linear 
ones in standard PT and correcting them by adding all possible 'chain-diagrams', i.e. 
diagrams in which 'soft' PS at momenta qi are attached to the hard line at momentum 
k, such that g, ^ k. The results above follow as a consequence of the property of the 
vertex function, eq. ([T]), 

1 k ■ q 

7a6c(k, -k + q, -q)Mc^5ab;T^-, for A; > g , (A. 3) 

Z q 

and of the composition property of linear propagator, 

9ab{s, s')gUs', s") = 9{s - s')9is' - s")gUs, s") . (A.4) 

Notice, in particular, as also discussed in the text, that the equal-time PS is equal to 
the linear one. 

These expressions are exact in the deep eikonal limit, that is for very large 
momentum k (neglecting multi streaming) and under the assumption that the 
cosmological scenario presents a clear hierarchy of scales allowing us to set apart 'soft' 
from 'hard' modes. Corrections to this idealized situation come from loop corrections, 
once they are properly taken into account. This can be done by using the path-integral 
formulation introduced in [25] (see also [22] )• In that paper, it was shown that standard 
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PT is equivalent to an expansion in powers of the vertex function 'jabc of the generating 
functional 

+ i J dl] [XaQab'^b - e'' 'JabcXa'fb^c + Ja^a + K^Xb] } , (A.5) 

where Ja and are sources for and Xb respectively, and we have omitted the 
momentum dependence. In deriving the above expression we have assumed Gaussian 
initial conditions. Non-Gaussian initial conditions can be taken into account by 



including a non-vanishing bispectrum, trispectrum, etc., in the first line of Eq. (A.5) 



Derivatives of Eq. (A.5) w.r.t. the sources Ja and Ki, give all the possible statistical 



correlators involving the fields and Xb-, such as the nonlinear PS and nonlinear 
propagator, given by 

1 5^Z 



Z<5J,(k,r7)5J,(k',r7') 
1 5^Z 



Ja,Kt = 



ZSJaiKT]) 6Kb{k',7]') 



SDik + k')Pab{k-r],r]'), 
-i 6D(k + k')Gab{k;r],T]') 



(A.6) 



Ja,Kt=0 

respectively. 

In order to pass from standard PT to eRPT, without modifying the full dynamical 
content of the generating functional, we simply add and subtract the quadratic 
expression 



dr]dr]' - lXa{v)Kb{v, v')Xb{v') - iXa{v)Kb{r], 



(A.7) 



to the exponent of eq. (A.5). Then, we include the added term in the new "free" {i.e. 



quadratic) part of the action, whereas the subtracted term goes in the new "interaction" 



part, which also includes the trilinear vertex at the third line of eq. (A.5). Of course 



since we have just added and subtracted the same term, standard PT and eRPT have 
the same dynamical content, and their results fully coincide at infinite loop order. 

The 'self-energy' and 'mode coupling' function in the eikonal limit, r/, r/') 

rj, rj'), can be represented by diagrams in which the hard line, carrying momentum 
of order k, is corrected by attaching soft PS in all possible ways such that the final 
diagram is one-particle-irreducible (IPI), that is, it cannot be cut into two disjoint pieces 
by cutting a single line. In terms of eikonal propagators, PS, and renormalized vertices, 
they are given by the 1-loop diagram of fig. [8| and by the 1- and 2-loop diagrams 
at the RHS of fig. [6} respectively, in which the tree rightmost vertex is replaced by 
a renormalized one (thick dot). In all these diagrams the propagator at the right, 
connecting s' to rj, has to be cut-off. 

Setting to zero the new interaction term, the path integral can be performed 
analytically, to give 

Zf^^[Ja, if,] =exp|- J dvdr^' [lUv)P:t{k;v,v)Jb{v') 
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propagator tree-level : -iG^ab ^^''''la,ilb) 



power spectrum tree-level: P^^^i'k;'qa,'rib) 



A 



interaction vertex: -ie^Yabc(^a,^b,^c) 




counter term propagator: il.^j-^(k;ria,rih) 



o 



a counter term power spectrum: -O^^^^Ck;?;^,?;^) 



Figure Al. The Feynman rules for eRPT. 



(A. 



where C^^^ and P^l^ are given in eqs. (A.l) and (A. 2), respectively, and are related to 
S^f and by 

= 9abiv,v') + J dsds'gac{yi,s)T.^^{k-s,s')Gt{k-s\r^')ik.Q) 

and 



(A.IO) 



The new expansion can be expressed diagrammatically as illustrated in fig. [Appendix A[ 
Besides the eikonal propagator and PS given above, in the interaction sector we have the 
usual trilinear interaction and two new counterterms, iS^^^f (fc; s, s') and — s, s'). 
These counterterms avoid overcounting of the contributions already included in the 
resummed propagator and PS. 

To see how eRPT works, we consider the propagator up to 1-loop order. The 



corrections to the tree level expression, (A.l), are given by the 1-loop contribution of 
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k . 

n' 




B 



Figure A2. The first non-trivial corrections to in eRPT 



fig. A2 



where 



^<lf(fc;s,s')7/5/.(|k-qU,fc) 



(A.ll) 



dg 

and by the contribution of the counterterm, 

= eal [ dse^e^ - e^')GZ'{k;v,s)Gll'{k;s,v') 



(A.12) 



v^^(e'-e''')Erf 



exp 



9ab{.ri, 7]') ■ 
(A.13) 



Notice that the 1-loop term, eq. (A.ll), contains the resummed propagator and PS, but 



tree level vertices. The latter appear in their exact form, that is, not in the eikonal limit 

(A. 14) 



of eq. (A. 3). In the small k limit the tree level propagator, eq. (|A.1|) goes as 



The 0{k'^) term is exactly canceled by the counter term in (A.13), which gives 



k'<y r ^' ' gabiv,v') + Oik') 



(A.15) 



(A.16) 



so that, summing up eqs. ( A.l[ ), (A.ll), and ( |A.13 ), one gets 

9abiv,v') + J ds ds'gac{v,s)J:^^,i{k;s,s')ghb{s',v') + Oik"^) , 
with 

El}^{k; s, s') = 4e^+^' J d'q ^^k, q, |k - q|)(7e/(s, s') x 

P%k)udUg-ffgh{\k-q\,q,k), (A.17) 
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that is, we obtain the propagator at 1-loop in standard PT, gabiv^^') + '^9abiv^''l')- 
order to get the 0(A;^) terms of standard PT one has to include the 2-loop diagrams 
of eRPT - not only the chain-like ones - and up to two insertions of the H'^j^ counter 
terms. Indeed, the single counter term insertion in (A. 13), expanded up to O gives 
—1/6 k^a^{e^ — e^'Y-. while the double insertion gives 

ds ds' ds" ds'" GZ\k; 7], s)Efa\k; s, s')G''l\k; s' , s") 

xi:if{k;s",s"')Gf{k;s"",v') 

gaM) + 0{k'), 



4! 



(A.18) 



so that the contribution from counterterms exactly cancels the k'^a^{e'^ — e''')^/8 term 
coming from the expansion of G'^^^. 

Notice that the appearance of the tree level vertex in its complete momentum 
dependence {i.e. not in the approximated form (A. 3)) is essential in order to recover 



the correct momentum dependence of the propagator, or, analogously, of the PS. 



The correction to the tree level propagator (A.l), i.e. the sum of (A. 11) and (A. 13) 



amounts to taking one of the soft PS insertion in the chain diagrams in the eikonal 
limit, and to replacing it with an insertion in which the vertices are at the tree level in 
their full form (instead of being fully renormalized, but in the eikonal approximation), 
and the inserted PS is also renormalized, i.e. it carries gaussian damping factors, as 



in (A.2). In other words, the contribution in (A. 11) represents a 1-loop "super-chain" 



diagram, in which, compared to the usual 1-loop chain diagram, all the propagators and 
the internal PS are renormalized. 

Differently from standard PT, at large k the different orders in this expansion do 
not diverge as some power law in k, but are always damped by the gaussian decay of the 
external eikonal propagators. In the large k limit the vertices can be approximated as in 



(A. 3) and the summation of the super-chain diagrams at all orders gives rise to a super- 



resummed propagator. In this resummation the counter terms should be properly taken 
into account in order to avoid over counting of the standard chain diagrams already 
included at the tree level. At intermediate scales, the full momentum dependence of the 
vertices, and also non-chain diagrams should be taken into account to properly compute 
the propagator away from the eikonal limit, as we have done up to 0(/c^) corrections. 
We now consider the PS. At 1-loop in eRPT we need to consider the diagrams in 



fig. A3, where also the counterterm now appears. Also in this case, the 1-loop PS 



in standard PT is reproduced, at low k up to O(fc^) terms. The role of the external 
propagator corrections, and of the corresponding counterterms T,*^^ is analogous to what 
we have just discussed for the propagator. As for the mode-coupling part, the sum of 



the tree level PS of eqs. (A.2) (see also (A. 10)) and of the counterterm gives the first 



line in eq. (A. 10), which is doubly gaussian suppressed, and plays therefore a sub leading 



role at intermediate scales. We are therefore left with the 1-loop (in eRPT) contribution 



dsds'GZHk;V,s)Glf{k;v' 



eik I 



(A. 19) 



Nonlinear Power Spectrum from Resummed Perturbation Theory: a Leap Beyond the BAO Scale33 



■Pi 

III 




k-q 



IV 



B k + k «i 

n n 




Figure A3. The first non-trivial corrections to P^l in eRPT 



where 



eRPT,{l) 
cd 



{k;s,s') 



2 e^+^' / d'ci^,,,{k, q, |k - q|)76e/(A;, q, |k - q|)P,f (g; s, s')Pf{\k - q|; s') , 

(A.20) 



Which we will use in eq. (49) in its approximated form valid for k ^ q. 



Appendix B. 



In this appendix we prove eq. (37). First, we notice that, as proved in [22], we have 
AG:f ' ^'\k-v,s)= ' ^'-'\k;r^,s) J' ds' J:Z''^'\k-v,s')u, , 

= H^ik;r],s)GZ'^^'~'\k;v,s), (B.l] 



where the loop indices refer to standard PT, not to the eRPT described in Appendix A 



Then, in order to reproduce (44), we have to show that a similar expression, but with 
the first factor replaced by ds' Tjae'^^\k; rj, s')ue, can be extracted from the third 
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Figure Bl. A generic contribution to A<I>^*j^'^' ^"'■'(fc; 77, s'). 



line in (37). Indeed, we can write 



eik I 



(n) 



The j'-th order contribution to the mode couphng function in the eikonal hmit 

" ad 



; 77, s') can be represented by diagrams in which the hne containing the hard 



hnear PS, carrying momentum of order A;, is corrected by attaching j soft PS in all 
possible ways such that the final diagram is one-particle- irreducible, that is, it cannot 
be cut into two disjoint pieces by cutting a single line. Then, we do the splitting 



eik, L, (j) 
ad 



{k;r],s') + 



eik, L, (j) 
ad 



{k;T],s') 



(B.3) 



where $ 



eik, L, (j) 
ad 



{k; 7], s') contains all the diagrams in which the soft PS attached to the 
left end, here identified by index "a" and time "77" , has the other leg attached to the right 

[k] rj, s') contains all the other diagrams. Equivalently, 



of the hard PS, and A$ 



eik, L, (j) 



the diagrams belonging to ^'^^^'^'^''\k;ri, s') have the leftmost vertex at the tre e level. 



Bl 



In the 



A generic diagram contributing to A^'^^^'^'^''\k;T], s') is given in fig. 
eikonal limit the soft PS attached to the leftmost end gives a factorized contribution of 
the form T,ae'^^\k; rj, s") Ug = —Ua k'^a'le^^'^" . This is a consequence of the property of 
the elementary vertex of eq. ( A.3 ), and of the composition property of linear propagator, 



eq (A.4) 



We are left with a contribution to the (j 



ds 



GZ'{k;v,s)^-'{k,s,s' 



eik I 



- l)-th order quantity 

(i-i) 



(B.4) 



where the time s corresponds to the point where the IPI function (^'^^^{k, s, s') is joined 
to the connected one, G^^{k;ri, s). Notice, in particular, that the time argument s" of 
the factorized S*^*^' ^^^-function is, in any case < s, since, otherwise, the corresponding 



diagram for A$ 



eik, L, (j) 
ad 



{k; 7], s') would not be IPI. Now, summing up all the contributions 



-^eik,L, (j)^ 

with the end of E^**'' at s" attached in all possible ways, gives 



to A$^^'*'''^' (A;; ?7, s') containing the same particular diagram contributing to eq. (B.4) 



ds / ds'EZ''^'\k-v,s')u, GZ'{k;v,s)<^lf{k,s,s' 



(B.5) 



Inserting (B.l) in the first line of ( [37| ) and ( B.5[ ) (and the corresponding contribution 
for a -H- b) in the third line, proves the equality. 
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Appendix C. 

In this appendix we will give the explicit expressions of the IPI functions Ti^^fj {k;ri,rj') 
and ^l^ij {k;r],r]') computed at 1-loop. The general expression for T^^^i^{k;r],r]') is given 
by 

4:6"^^"' j d^q-facd{K -q, q - k)McP°(g)Me7/efe(k - q, q, -k)^d/(r/, VlC^.l) 

which, using the expressions for the vertices and the linear propagator in eqs. ([T]), (14), 
and performing the angular integration, gives 



-3/2r?+r; 



~15 



dr P^kr) 



+ e^/^^\l7 + 9r2) + |- (e^/^'? _ e^/^'?') (r^ _ i)2 log ^ + 



E«(fc;r/,V) = -e-'/'''+^'^(e^/'^-e^/2''') J dr P%ky 



1 — r 

.2 I ri»,4 



8 - ISr^ + 9r^ 



1 +r 



1 — r 



s(V(A;;r/,r7') = _e-3/2^+V^^(e5/2r, _ ^5/2,') /" ^^po^^^) rj_^3^2 _ 

5 J [r^ 



1 + r 



1 — r 

sg(A;;r/,V) = j dr P\kr) 



(C.2) 



1 


_l_ ^ 




1 


— r 





_(e5/2._e5/2,'))|_^^2_^^2iQg 

The general expression for the 1-loop contribution to ^^^^{k^rj^rj') is 
$11)(A;;77,V) = 

2e''+^' j d^q^acdiK -q, -p)ucP°(g)neUdP°(p)u/7;,e/(-k, q, p) , 
with p = k — q, which gives the following expressions for the individual components. 



(C.3) 



Ak 



Jo J\k-q\ P 1 



4 Jq J\k-q\ P-^T 
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\k-q\ 



p3q3 



-P\q)P\p) 



(C.4) 



Appendix D. 

In this appendix we will give the explicit formulae to be inserted in the evolution equation 
for the nonlinear PS, eq. (53). The first line contains the matrix Vtab-, defined in eq. ([9]). 
At the second line, the two functions Hi{k; 77, —00) and H2{k; rj, —00) appear, which are 
given explicitly by 

Hi{k;r],-oo) = / ds {k; rj , s)ub = 



— e 



— I dr 
21 



19 - 24r2 + 9r^ - 



(r^ - if lo£ 



H2{k;r], -00) 



2r 

dsJ:''^^{k;ri,s)ub 



1 + r 



1 — r 



P%kr) 



(D.l) 



~2l 



dr 



2r3 



+ 52 + 9r^ 



l)'log 



1 + r 



1 — r 



P^kr) 



(D.2) 



As we have discussed in the text, the third line of eq. (53) is well approximated by 



using the small scale expression for given in eq. (51), and the small scale analytic 



expression for the propagator, G^^, instead of the solution of eq. (25) with eq. (28) 
namely. Gab- The time integration can then be done analytically, to get. 



ds [^adik; 7], s)G^^{k; t], s) + G^^ik; r], s)^dbik; s, 7])] 

= ^Gi{k;r^) + ^G^,{k;v), 



(D.3) 



where 



^Gt,{k;v) 



7rErf(y) ( 3$SlH^;^,^i ^ 2^]"^ {k; 7] , r]) 



y \ 5 ' 5 

+ ^B{y'){4\\k;v,v)-^^S{k;V,v)) 



^GUk;v) 



¥Erf(y) /3$S(/^;^,^) , '^\'^ik;V,v) , ^'22ik;v,v) 



y 



10 



+ 



+ 



-—Biy') (3<l>[V(fc;r7,r/) - 5<l>f2^(fc; r^, r/) + 2<|.g(A;; r/, r/)) 
5/5 

^G^,{k;r]) = ^GUk;v), 

22 ' y 
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Figure Dl. The results in the BAO region for k = 0.15, 0.20, 0.25, 0.30 hMpc"\ 
(orange, red, cyan, purple lines, respectively) 
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and 



^G'!,{k-7^)=UaU,P\k) 



1) 



2/'Erf(y) 



(D.4) 



(D.5) 



with y = e^kcTy. The function B{y'^) is a combination of generahzed hypergeometric 
functions, 



W =35 2^2 (-,1;-, 



1 L 

'4'4' 
3 . 9 11 



28 2F2(i,|^,^;-r 



+ 82^2! ^,2;-,-; -r 



(D.6) 
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Eq. ( D.3[ ) gives the large y limit of t he t hird line of e g. (|53| ). In the very large k limit 
^^ab i^'iV^v) ~^ UaUbP^{k)y'^ and eq. (D.3) goes to eq. (55). 

As we have discussed in the text, in the small momentum limit the second term, 
$Gf^(/c; rj) has to be switched off, because it contains 2-loop expressions valid at large 
k. Therefore we will multiply it by a momentum cutoff function. In our numerical 
implementations, as discussed in the text, we have used the expression 

,4 



{k;r]) 



(I) 



(D.7) 



where k = 0.2h/Mpc represents a reasonable value above which the large scale 
expression can start to be applicable. Alternatively, k can also be taken as a parameter 
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to be marginalized in fits to real data. In figs. |D1| and |D2| we show results obtained by 
changing the filter value from k = 0.15 to 0.30 hMpc~^. 

The initial conditions must be given at a very large redshift Zm = O(IOO), where 
the PS can be approximated with the linear one: 

Pai>ik;r],^) = UaU,P''ik), (D.8) 

and the equations can then be integrated town to the required final redshift, 
corresponding to rjf = log D{z f) / D{zin). 
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